DISSIPATIVE PROPERTY OF THE VLASOV-MAXWELL-BOLTZMANN 
SYSTEM WITH A UNIFORM IONIC BACKGROUND 

RENJUN DUAN 

o 

^N) ' Abstract. In this paper we discuss the dissipative property of near-equilibrium classical solu- 

tions to the Cauchy problem of the Vlasov-Maxwell-Boltzmann System in the whole space R^ 
when the positive charged ion flow provides a spatially uniform background. The most key point 
of studying this coupled degenerately dissipative system here is to establish the dissipation of 
the electromagnetic field which turns out to be of the regularity-loss type. Precisely, for the 
linearized non-homogeneous system, some L^ energy functionals and L^ time-frequency func- 
tional which are equivalent with the naturally existing ones are designed to capture the optimal 
dissipation rate of the system, which in turn yields the optimal U'-L'^ type time-decay estimates 
P^ ' of the corresponding linearized solution operator. These results show a special feature of the 

^^ ' one-species Vlasov-Maxwell-Boltzmann system different from the case of two-species, that is, the 

^J dissipation of the magnetic field in one-species is strictly weaker than the one in two-species. As 

(*^ , a by-product, the global existence of solutions to the nonlinear Cauchy problem is also proved by 

constructing some similar energy functionals but the time-decay rates of the obtained solution 
still remain open. 
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1. Introduction 

1.1. Main results. The Vlasov-Maxwell-Boltzmann system is an important model for plasma 
physics to describe the time evolution of dilute charged particles (e.g., electrons and ions in the case 
of two-species) under the influence of the self-consistent internally generated Lorentz forces [TB] . In 
physical situations the ion mass is usually much larger than the electron mass so that the electrons 
move much faster than the ions. Thus, the ions are often described by a fixed ion background 
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nb{x) and only the electrons move. For such simple case, the Vlasov-Maxwell-Boltzmann system 
takes the form of 

f aj + C • V,/ + (i? + e X B) • Vj/ = Q{f, /), 



(1.1) 



dtE-V^x B = - 
dtB + W^xE^ 0, 



V -E 



fd^-n^, V,-B^Q. 



Here, the unknowns are / = /(i, a;, ^) : (0, oo) 



[0,oo), E :^E{t,x) : (0,oo)x 



and B = B{t,x) : (0, oo) x R'^ -> R"^, with f{t,x,^) standing for the number distribution function 
of one-species of particles (e.g., electrons) which have position x = {xi,X2,X3) and velocity ^ = 
(CijCzjCs) at time t, and E{t,x) and B{t,x) denoting the electromagnetic field in terms of the 
time-space variable (i, x). The initial data of the system at i = is given by 

(1.2) /(O,x,0 = /o(a:,0, EiO,x) ^ E^ix), B{0,x) ^ Bo{x). 

Q is the bilinear Boltzmann collision operator [1] for the hard-sphere model defined by 



Q(/,.9) 



'XS2 



{f'g* - ./.9*)I(C - ■?*) • uj\dujd£,^, 



f^f{t,x,0, f'^fit,x,C), g.^9{t,x,^,), 9:^9it,x,C), 

^' = c-[(e-e*)-Hw, e: = e* + [(e-c*)-w]w, coes^ 

Notice that system (jl.ip in general contains physical constants such as the charge and mass of 
electrons and the speed of light, cf. [11]. Since our purpose in this paper is to investigate the 
dissipative property of solutions near global Maxwellians, those physical constants in system (jl.ip 
have been normalized to be one for notational simplicity. Through this paper, ni,{x) is assumed 
to be a positive constant denoting the spatially uniform density of the ionic background, and we 
also set TT-b = 1 without loss of generality. 

We are interested in the solution to the Cauchy problem for the case when the number dis- 
tribution function f{t,x,£^) is near an equilibrium state M and E{t,x) and B{t,x) have small 
amplitudes, where M denotes the normalized Maxwellian 

M = M(0 = (27r)-3/2e-l?lV2. 
For that, set the perturbation u by 

/(i,x,0 = M + Mi/2^(t,x,e). 
Then, the Cauchy problem (jl.ip . (|1.2p can be reformulated as 



(1.3) 



dtu- 

dtE 
dtB 



- C • V^u + {E + ^x B)- W^u - CM^/2 • E 

= Jju + r(u, u) 



-f ■ Eu, 
2^ 



Vj; X S 

V^E X £■ 



SM^/^udS,, 



\/ ■E= [ M^/^ud^, S/x-B 

JR3 
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with initial data 

(1.4) uiO,x,0 = Mx,0, EiO,x)=Eoix), BiO,x)^Boix). 

Here, the hnear term Lu and the nonhncar term r{u,u) are defined in (|2.1[) and (|2.2p later on. 
The problems to be considered are (i) whether or not any small amplitude solution 

[u{t),E{t),B{t)] ■.R+ ^ X = H^iRl X M|) x -ff^(M^) x H^ (Rl) 

with a properly large N for the above reformulated Cauchy problem uniquely exists for all t > if 
initial data [uq, Eq, Bq] € X is sufficiently small; (ii) if so, does the solution decay in time with some 
explicit rate? We shall give in this paper a satisfactory answer to the first question and a partial 
answer to the second one only in the linearized level. Since these two issues have been extensively 
studied in different contents such as the Boltzmann equation [H [3] , the Vlasov-Poisson-Boltzmann 
system [71 [B] and the Vlasov-M axwell-Boltzmann system for two-species [S] , our emphasize here 
will be put on the study of the weakly dissipative property of the electromagnetic field and its 
resulting slow time-decay rate of solutions, which are even different from the case of two-species. 
All details for their discussions are left to the next subsection. 

Let us begin with the Cauchy problem on the linearized non-homogeneous Vlasov-Maxwell- 
Boltzmann system, in the form of 



(1.5) 



dtE^W^x B^- ^M^/^ud^, 
dtB + V,xE = 0, 

with 

(1.6) u(0,x,0 = "o(a;,e), E{Q,x)=E^{x), B{{),x) ^ Bq{x), 

where h = h(t, x, ^) denotes a given non-homogeneous source term. For simplicity, we write 

U=[u,E,B], Uo^[uo,Eo,Bo]. 

Moreover, Uq — [uq, Eq, Bq] is always supposed to satisfy the last equation of (|1.5p for t = 0. The 
first result, concerning the naturally existing energy functional and its optimal dissipation rate, is 
stated as follows. Here and hereafter, i/ = i/{^) and P arc defined in p.3|) and ()2.5p . respectively; 
see Subsection 11.31 for more notations used in this paper. 

Theorem 1.1. Let N > 3. Assume v^^^^h £ Lj{H^) with Ph{t,x,£,) = 0. Define the temporal 
L'^ -energy functional £^^^{U{t)) and its dissipation rate T)^™ [U {t)) by 

(1.7) £'^{U{t)) ^ lh^(t)!li.(^«) + ||[i?(i),i?(t)]||^„, 

(1.8) v's'iuit)) = \\u^r2{i _ v}u{t)\\l.^^^^^ + ||v.Pu(i)|l2,^^„_,^ 



2 I IIV72 D,'j.M|2 



\y,E{t)r^r.-. + \\^iB,.,,^ 



Ar-3- 



Then, for any smooth solution U ~ [u,E,B] of the Cauchy problem (|1.5p - (|1.6p belonging to 
L'j{H^) X H^ X H^ , there exists a continuous functional £^™ {U {t)) given in p.l2p such that 

(1.9) 1^^"(^(^)) + APl^"([/(t)) < C\\u-'"h{t)\\lii^HS) 

for any t > 0. 
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Remark 1.1. The above theorem shows the precise dissipative property of the naturally existing 
L^ -energy functional 

\\uit)\\l.^H^^ + \\[Eit),B{t)]\\l^ 

for a properly large N. The construction of the equivalent energy functional £]^"(C/(i)) is used to 
capture the optimal dissipation rate T)^™{U{t)), which will also he revisited in Corollary \1.2[ This 
is different from the case of the two-species Vlasov- Maxwell- Boltzmann system as in [S|, where the 
dissipation rate is a little stronger due to a cancelation phenomenon between two species which 
was firstly observed in [T7]. Specifically, if it is here supposed that h = and E]^{Uo) is finite, 
then not only all macroscopic quantities Pu, E, B and the highest-order derivative V^[£',-B] of 
the electromagnetic field lose their time-space integrability , but also the same thing happens to the 
second-order derivative V^i? of the magnetic field. On the other hand, for the energy space with 
m-order spatial regularity for any integer m > 0, its optimal dissipation rate can be described by 
I?Jj"([/(i)) once again from Corollary \1.2[ Finally, as seen from Theorem \1.S\ later on, this kind 
of weaker dissipation property leads to some slower time-decay rates of solutions. 

The second result about some time-frequency functional and its optimal dissipation rate is stated 
as follows. 

Theorem 1.2. Assume v^^^'^h e L? for t > 0, k £ M.^ , and Ph = 0. Define the L^ time-frequency 
functional £{U{t, k)) and its dissipation rate T>{U{t, k)) by 



(1.10) £{U{t,k))^ Ml, + \[E,B]\\ 



(1.11) v{u{t, k)) = ||z.V2{i _ p}^||^, + |fc . ^|2 + ^i-L_|[a, b, 

\k? ,^,2 . \k\' ,5,2 



c]? 



EV + ,_ ' ,,..J B 



(l+|fc|2)2' ' (l + |fc|2)3 

Then, for any solution U — [u, E, B] of the Cauchy problem (|1.5|) and (|1.6p satisfying that \\u\ 



2 



\[E, i?]|2 is finite for i > and k G R'^, there indeed exists a time-frequency functional £{U(t, k)) 
given in (|3.30|) such that 

(1.12) dt£{U{t, k)) + XV{U{t, k)) < C\\v-^/^h\\l2 

for any t > and fc G M"^. 

Corollary 1.1. Under Theorem \1.2\ £{U{t,k)) further satisfies 

(1.13) dt£{U{t, k)) + ^-Mj_£(f7(i, k)) < CJk-i/2/,||2 , 

for any t > and fc G M'^. 

Remark 1.2. Theorem \1.2\ and Corollary \1.1\ show that the magnetic field B bears the weakest 
dissipation property among all quantities {I — P}u, a, b, c, E and B, and even the dissipation of 
B here is much weaker than that in the case of two-species Vlasov-Maxwell-Boltzmann system as 
in [5]. This is also consistent with what has been mentioned in Remark ] 1. 11 
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Corollary 1.2. Let m > be an integer. Assume v^^l'^\I'"^h G i^ with Ph{t,x,S,) = 0. Define 
the L^ energy functional E]™ [U [t)) and its dissipation rate T)^™ {U {t)) by 

(1.14) f;;r(^W)-llv>(i)lP + l|vr[i?W,i?W]|p, 

(1.15) vll^{u{t)) = \\^'/^w":{i-p}u{t)r 

+ l|V™a|r + ||Vi+™(V.)-Ma,&,c]f 

Then, for any smooth solution U ~ [u, E, B] of the Cauchy problem (|1.5p and (jl.6p whose m-order 
spatial derivative belongs to L^ ^x L^ x L^ , there indeed exists a continuous functional 8^^{U{t)) 
given in p.3ip such that 

(1.16) 1^™ (^W) + ^-D'^^Uit)) < C\W-'/'V^hit)r 

for any t > 0. 

Remark 1.3. It is straightforward to observe that Corollary \1.S\ imvlies Theorem \l.l\ by defining 

N N 

^N\U{t)) = J2 ^™ (C^W), T^NiU{t)) - J2 ^™ (f/w), 

m— m— 

and using 

N N 

E iv.r+'"(v.)-i = iv.i(v.)-i Y. i^-r - iv.i(v.)^-s 

rn—O rn—0 

and likewise 

N N 



J2 |Vxr+™(v,)-2 ^ |v,|(v,)^-2^ 5] |V,|2+"(V,)-3 ^ |V,|2(V,)^ 



Notice that the above identities and equivalent relations can be verified with respect to the frequency 
variable under the Fourier transform. On the other hand, it is also interesting to see that even when 
< m < N, I?Jjj"(t/(t)) can capture the optimal dissipation rate of the naturally existing m-order 
energy functional £l^{U{t)) . For instance, when m ~0, the direct energy estimate on system (|1.5p 
produces the only dissipation of the microscopic component {I — P}u, which is partially contained 
in the optimal form 2?g"([/(t)). 

Furthermore, we can obtain the large-time behavior of solutions to the linearized non- homogeneous 
Cauchy problem. Formally, the solution to the Cauchy problem (jl.Sp and (|1.6p is denoted by the 
summation of two parts, 

(1.17) U(t) = U'{t) + U"{t), 

(1.18) U' (t) ^ A{t)Uo, U' = [u^,E',B'], 

(1.19) U"{t)= f A{t-s)[h{s),0,0]ds, U^^ = [u",E^^,B^^], 

JQ 

where A(i) is the linear solution operator for the Cauchy problem on the linearized homogeneous 
system corresponding to (jl.Sp with h — Q. Notice that U^^{t) is well-defined because [h{s),0, 0] for 
any < s < i satisfies the last equation of (jl.Sp due to the fact that Ph{s) — and hence 

/ M^/^h{s)d£.=0. 
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For brevity, we introduce the norms || • \\h"^, \\ • Wz^ with ttt, > and r > 1 given by 

for U = [u, £', B], and note Z2 = Ti.^ ■ Then, the third resuh to describe the time-decay property 
of the hnearized solution is stated as foUows. 

Theorem 1.3. Let I < p,r < 2 < q < 00, cr > 0, and let m > be an integer. Suppose Fh = 0. 
Let U he defined in (|1.17p . (|1.18p and (|1.19p as the solution to the Cauchy problem (jl.5p - (|1.6p . 
Then, the first part U^ corresponding to the solution of the linearized homogeneous system satisfies 

(1.20) ||V™[/^(i)b, <C(l+i)-5(i-f)-T||f/o||^^+C(l + r^l|Vr'''''''^"'"+t/o|U,, 

for any t > 0, and the second part U^^ corresponding to the solution of the linearized nonhomoge- 
neous system with vanishing initial data satisfies 

(1.21) ||v;"c/^^(i)|||, <c Ai + f-.)-i(^-^)-™||^-i/2Ms)lll,d^ 

Jo 

+ C [ {I + t - s)-^\\v-^'''^Z'+^''^+h{s)f ds, 
Jo 

for any t > 0. Here, [•]+ is defined by 

11 I (7 if a is integer and r = q = 2, 

(1.22) [a + 3( )]+ = 



^ I [o- + 3(^ - -)] + 1 otherwise, 



where [■] means the integer part of the nonnegative argument. 

Finahy, let us go back to the Cauchy problem on the nonlinear Vlasov-Maxwell-Boltzmann 
system. The global existence and uniqueness of solutions are stated as follows. 

Theorem 1.4. Let N > 4. Define L^ energy functional En {U {t)) and its dissipation rate'DN{U{t)) 
by 

(1.23) £N{U{t)) ^ \Ht)\\l. + \\[E{t),B{t)]\\]j., 

(1.24) VN{U{t)) = ||z.i/2{i _ P}u(i)||^„ + ||i.i/2v^y(t)||2 

x,^ ^^\^x ) 

+ ||V.i?(t)||^„_.+||V^i?(t)f^„_3. 

Suppose /o = M + M-'^'^Mo > 0. There indeed exists a continuous functional £]\[{U{t)) given in 
(j4.7p or (j4.8p such that if initial data Uq = [uq, Eq, Bq] satisfies (jl.3P zi for t = and £m{Uo) 
is sufficiently small, then the nonlinear Cauchy problem (jl.3p - p.4p admits a global solution U = 
[u, E, B] satisfying 

f{t,x,0 = M. + m^'^u{t,x,i)>Q, 

[u[t),E{t),B{t)] e C([0, 00); iJi^^ X H^ X ijf ), 

and 

£NiU{t)) + \ f VN{U{s))ds < £n{Uo) 
Jo 
for any t > 0. 

The decay rate of the solution obtained in Theorem 11.41 remains open. We shall discuss at the 
end of this paper the main difficulty of extending the linear decay property in Theorem 1 1.31 to the 
time-decay of the nonlinear system. 
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1.2. Related work and key points in the proof. As mentioned before, the method of con- 
structing energy functionals or time-frequency functionals to deal with the global existence and 
time-decay estimates presented in this paper has been also extensively applied in [H [71 HI 121 [5] . 
Specifically, jj is a starting point of these series of work. In U, some interactive energy func- 
tionals were constructed to consider the dissipation of the macroscopic part of the solution and 
also the global existence of solutions without any initial layer was proved. Later, the same thing 
was done in [7] for the one-species Vlasov-Poisson-Boltzmann system, where the additional ef- 
forts are made to take care of the coupling effect from the self-consistent potential force through 
the Poisson equation. In order to investigate the optimal rate of convergence for the one-species 
Vlasov-Poisson-Boltzmann system, a new method on the basis of the linearized Fourier analysis 
was developed in [6] to study the time-decay property of the linear solution operator, where the 
key point is again to construct some proper time- frequency functionals so as to capture the optimal 
dissipation rate of the system. At the same time, |3j provided another method to study the expo- 
nential time-decay for the linear Boltzmann equation with a confining force by using the operator 
calculations instead of the Fourier analysis. 

Recently, following a combination of [^ and |17j . the optimal large-time behavior of the two- 
species Vlasov-Maxwcll-Boltzmann system was analyzed in [5]. The main finding in [5] is that 
although the non-homogeneous Maxwell system conserves the energy of the electromagnetic field, 
the coupling of the Boltzmann equation with the Maxwell system can generate some weak dissipa- 
tion of the electromagnetic field which is actually of the regularity-loss type. It should be pointed 
out that even though the form of two-species Vlasov-Maxwell-Boltzmann system looks more com- 
plicated than that of the case of one-species, the study of global existence and time-decay rate 
is much more delicate in the case of one-species because the coupling term in the source of the 
Maxwell system 



JS.3 



corresponds to the momentum component of the macroscopic part of the solution which is degener- 
ate with respect to the linearized operator L. Essentially, it is this kind of the macroscopic coupling 
feature that leads to some different dissipation properties between two-species and one-species for 
the Vlasov-Maxwell-Boltzmann system. 

For the convenience of readers, let us list a table below to present in a clear way similarity 
and difference of the dissipative and time-decay properties for three models: Boltzmann equation 
(BE), Vlasov-Poisson-Boltzmann system (VPB) and Vlasov-Maxwcll-Boltzmann system (VMB). 
In Table 1, 1-s means one-species and 2-s two-species. Corresponding to different models, £(t,k) 
stands for some time-frequency functional equivalent with the naturally existing one and 2?(i, k) 
denotes the optimal dissipation rate of £{t, k) satisfying 

—£(t,k) + XV(t,k) <0 
at 

for all i > and fc G M^. All estimates are written for the linearized homogeneous equation or 
system. For more details and proof of other models in the above table, interested readers can refer 
to [3l [6l [5] . Notice that the 1-s VMB system decays faster than the 1-s VPB system due to the 
choice of initial data, that is, Eq £ L]. is assumed in the case of 1-s VMB system, whereas in 
the case of 1-s VPB system, the electric field Eq = —'^x'Po with the potential force (/)o satisfying 
the Poisson equation —AxCJ^o = /k3 M^/^mq d^ ni^Y "^ot belong to L], under the assumption of 
uq G Zif] L^. We remark that the decay rate (1 -I- i)~^/^ for the 1-s VPB system can be improved 
to be (1 + t)^'^/"^ provided that Va;(/'o G Ll. is additionally supposed. Finally, it should also be 
pointed out that the method developed in O [6l [5] and this paper could provide a good tool to 
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deal with similar studies for other physical models with the structure involving not only the free 
transport operator but also the degenerately dissipative operator; see also [20] , 



£{t,k) 



V{t,k) 



Ht)\\L^ < 



BE 



r 2 



^l/2{I_p}^||2 



C{l + t)-i\\uo\\z,n 



L2 



1-s VPB 



W 



.l/2|I_p}^||2 



"l+|feP 



[a,b,c]\^ + \a\' 



Cil + t)—^\\uo\\z,n 



L2 



1-s VMB 



i^ 



'/Hi-p}uh. 



\[E,B]\^^ 






c(i + t)-t 

x(||;7obi + ||V,[/o|U, 



2-s VMB 



\[E,B]\' 



z.V2{i_p}^|| 



[a±,b,c] 



i+|fep 
"(i+|feP) 



ifc-£;i' 



^(I^P + l^n 



c(i + t)-t 



Table 1 . Dissipative and time-decay properties of different models 



Since the current work is a further development in the study of the Vlasov-MaxwcU-Boltzmann 
system as in [5], we omit the detailed literature review for brevity, and readers can refer to [5] 
and reference therein. Here, we only mention some of them. The spectral analysis and global 
existence for the Boltzmann equation with near-equilibrium initial data was given by [19j . For 
the same topic, thirteen moments method and global existence was found by [13j . The energy 
method of the Boltzmann equation was developed independently in ^ [TUl E] and |15l UM \n\ 
by using the different macro-micro decomposition. The almost exponential rate of convergence of 
the Boltzmann equation on torus for large initial data was obtained in [2] under some additional 
regularity assumption. |18j provided a very simple proof of [2] in the framework of small pertur- 
bation. The diffusive limit of the two-species Vlasov-Maxwell-Boltzmann system over the torus 
was studied in [5j. It could be interesting to consider the same issue as in [S] for the one-species 
Vlasov-MaxwcU-Boltzmann system because of its weaker dissipation property. 

In what follows, let us explain some new technical points in the proof of our main results which 
arc different from previous work. The first one is about the dissipation estimate on the momentum 
component b{t, x) in the macroscopic part Pu in Thcorcm ll.il Recall that in the previous work [7] 
and [3], the dissipation estimate of b{t, x) was based on ()2.7P i and (|2.6|) 9. This fails in the case of 
one-species Vlasov-Maxwell-Boltzmann system because it is impossible to control a term such as 



E' 



c{diEj+djEi)dx. 



Instead, the right way is to make estimates on ()2.9|) zi. Therefore, one can get the macroscopic 
dissipation estimate p.8|) with the dissipation of E multiplied by a small constant on the right- 
hand side. The second key point is about the dissipation of the electromagnetic field. Different 
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from |5], although there is no cancelation in the case of one-species, one can still design some 



interactive energy functional £^' {U{t)) to capture the weaker dissipation rate 



l<|Q|<Af-l 2<|a|<iV-l 

The third key point is about the time-decay estimate on the linearized solution operator. In 
fact, for a general frequency function 4>{k) which is of the regularity-loss type as in Lemma 13.11 
one can repeatedly apply the Minkowski inequality to make interchanges between frequency and 
space variables so that the more general U'-L'' type time-decay than in Theorem 11.31 can be 
obtained. The last key point is about the estimate on the nonlinear term ^ • Eu in the proof 
of Lemma 14.11 concerning the a priori estimates of solutions. In particular, it is impossible to 
bound //^ • EiV'uf' dxd^ by using EN{U{t)Y/'^Vis[{U{t)) up to a constant since both E and Pm 
do not enter into the dissipation rate T>N{U{t)) given in (|1.24p . Instead, we first take the velocity 
integration and then use the macroscopic balance laws ()2.6p so as to obtain an estimate as 



£. ■ E(Pu? dxd£_ < — \b\^{a + 2c)dx + C £NiU(t))^^^ + £N{U(t)) VN{U{t)). 
dt J L J 

We remark that this also has been observed in [7| in the study of the one-species Vlasov-Poisson- 
Boltzmann system. 

1.3. Notations. Throughout this paper, C denotes some positive (generally large) constant and 
A denotes some positive (generally small) constant, where both C and A may take different values 
in different places. In addition, A ^ B means \A < B < jA for a generic constant < A < 1. 
For any integer m > 0, we use -ff"V, H"^, HV^ to denote the usual Hilbert spaces iJ'"(Ri^ x R5), 
i/™(M^.), iJ™(K5), respectively, and L^, L^, L? are used for the case when m = 0. When without 
confusion, we use if™ to denote H™ and use L^ to denote Li^ or L| t. For a Banach space X, |j • \\x 
denotes the corresponding norm, while || • || always denotes the norm || • \\]^2 for simplicity. For r > 1, 
we also define the standard time-space mixed Lebesgue space Zr = L'i{L^) = L^(M5; L'''(M.^)) with 
the norm 

\\9\\z,.= l[ ([ \9{x,0rdx] dA , g = g{x,OeZr. 

For multi- indices a = [ai,a2,a3] and /3 = [/3i,/32,/33], we denote 

The length of a is |a| — ai + a2 + 0.3 and the length of /3 is |/3| ~ /3i + (32 + Ps- For simplicity, we 
also use dj to denote dx for each j — 1,2,3. For an intcgrable function g : K'^ ^ M, its Fourier 
transform is defined by 

3 

where i = \/— 1 G C is the imaginary unit. For two complex vectors a, fe G C'^, (a | 6) denotes the 
dot product of a with the complex conjugate of b over the complex field. (V^) = (1 + |Vj;p)^/^ is 
defined in terms of the Fourier transform. 

The rest of this paper is arranged as follows. In Section [5] we present some basic property 
of the linearized collision operator and derive some macroscopic moment equations. In Section 
|3] we study the linearized non-homogeneous Vlasov-Maxwell-Boltzmann system in order to prove 
Theorem 11.11 Theorem 11.21 and Theorem 11.31 Finally, we prove in Section |4] Theorem 11.41 for the 
global existence of solutions to the nonlinear Cauchy problem. 



g[k)=Tg{k)= / e ^"^'''g{x)dx, x ■ k :=y^^j^j^ ^^ 



10 



RENJUN DUAN 



2. Moment equations 
It is easy to see that Lw and T{u,u) are given by 



(2.1) 
(2.2) 



LiU = 



1 



T{u,u) 



M 
1 



Q(M, VMu) + QiVMu, M) 
Q{VMu,VMu). 



For the hnearized cohision operator L, one has the fohowing standard facts [I]. L can be spht as 
JjU = —iy{S,)u + Ku, where the coUision frequency is given by 

(2.3) H0= ff |(e-e*)-^|M(^*)dc^dC,. 

Notice that i'(^) ~ (1 + 1^^)"^^^. The null space of L is given by 

AA = span|Mi/2,<e,Mi/2 (1 < i < 3),\^\^M^/^\ . 

The linearized collision operator L is non-positive and further — L is known to be locally coercive 
in the sense that there is a constant Aq > such that [T] : 

(2.4) - / uLudC>Xo I i'(0|{I-P}Mpde, 

where, for fixed (t,x), P denotes the orthogonal projection from L? to J\f. Given any u(t,x,^), 
one can write P in (|2.4I) as 

(2.5) Pu = {a{t, x) + b{t, x)-^ + c{t, x){\^\'^ - 3)}M^/^. 

Since P is a projection, the coefficient functions a{t,x), h{t,x) = [&i(t, a;), 62(^,2;), &3(i, a;)] and 
c(t, x) depend on u{t, x, ^) in terms of 

'a= / m^/'^ud£,^ I m^''^V>ud£„ 

Jm3 Jm.3 

< b,= [ CM^^^udC = [ ^M^^^FudC, 1 < i < 3, 
c^l f {\e-i)M'^^ud^=l f {\e-3)M'^'Pudt 

*- 6 Jr3 6 Jr3 

To derive evolution equations of a, b and c, we start from the local balance laws of the original 
system (|l.ip to obtain 



dt / /rfc + V, • / e/rfe = 0, 

JR3 JR3 JR3 JR3 

'^ JR3 ^ JR3 ^ JR3 



The above system implies 

dta + \7^-b = 0, 



(2.6) 



dtb + \/^{a + 2c) + \/^- I ^ (^ ^M^/2{j „ pj^f^^ - S(l + a) - 6 x B = 0, 

JR3 

a^c + iv, • & + iv, • / |^|2^Mi/2{I - P}udx -1^.6 = 0. 
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(2.7) 



As in [6], define 

0yH= / m^ - l)M'/'ud^, A,(u) = ^ / {\^f - 5)^M'/'ud^ 

JR3 lU Jr3 

for 1 < i,j < 3. Applying them to the first equation of (jl.3p . one lias 

' at[e.y({i - p}u) + 2cS,j] + d,b, + djh -= e,, (^ + g), 

where 5ij means Kroncckcr delta, and 

(2-8) <^ 1 

\g = -i-Eu-{E + ixB)- Vju + r(u, u). 

One can replace dtc in (|2.7l) i by using (HiUla so that 

2 

3 

10 



^tQ^J{{l - P}u) + d,bj + djh - :j<5,jV, • b 



- y%V, • A({I - P}u) = e,,(£ + ff) - -5,,-B • 6. 



In a summary, we obtained the following moment system 

dtb + \7^{a + 2c) + V^^ • e({I - P}u) -E = Ea + bx B, 



(2.9) 



dtC + i V, • 6 + ^V, • A({I - P}u) ^^E-b, 

2 10 

dte,,{{I - P}u) + d,bj + d,b, - g-JyV, . 6 - y%V, • A({I - V}u) 

2 
3' 



= e,; 



-g)^-S,jE-b, l<tJ<S, 



[ dtM{I - P}u) + d,c = K{i + g), 1 < z < 3. 

On the other hand, the Maxwell system in (|1.4p is equivalent with 

dtE^W.,xB^ -b, 

dtB + S/a;X E ^ 0, 

V, -^ = 0, V, -5 = 0. 

Finally, we should point out that the key analysis of all results in this paper is based on the above 
moment equations coupled with the Maxwell system. 

3. Linear non-homogeneous system 

In this section we consider the Cauchy problem (|1.5p and (|1.6p on the linearized Vlasov-Maxwell- 
Boltzmann system. For convenience of readers, recall it by 



(3.1) 



dtE~V^x B = - / CM^/'^ud^, 
dtB + V^xE = 0, 

V ■E= [ M^/^ud^, V^ • B = 0, 

JR3 
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with 

(3.2) u{0,x,O = Mx,O. E{0,x)=Eo{x), B{0,x) ^ Bo{x). 

Here, h ~ h{t, x, ^) is a given non-homogenous source term, satisfying Ph — 0. 

3.1. L^ energy functional and its optimal dissipation rate. In this subsection we shall 
prove Theorem ll.il Before that, similar to obtain (|2.9p . one can also derive the following moment 
equations corresponding to the linear equation p.l^ : 

r a^a + V, • 6 = 0, 

dtb + V^(a + 2c) + V^ • e({I - P}u) - E = 0, 



(3.3) 



dtc + iv., • 6 + ^V, • A({I - P}u) = 0, 

2 10 

9*6^ ({I - P}u) + d,hj + djh - 3<^,j V, • h - y%V, • A({I - P}u) = e,,(^ + h), 



ydtH{l-P]u)+d,c = k,{l + h), 
where 1 < «, j < 3, and as in (|2.8p . ^ still denotes 

£= -^. V:r -{I-Plu + Lu. 
The Maxwell system also takes the form of 

{dfE-W^x B ^ -b, 
dtB + V,xE^ 0, 
V, •£; = 0, V, -5 = 0. 

Proof of Theorem ll.lt Let A^ > 3. First of all, a usual energy estimate on (|3.ip gives 

(3-5) I E |(ii9%r + !i5"[i?,s]f)+A 5^ ii^v2aa{i_p}^||2<c ^ 11^-1/25, 

|Q|<Af |Q|<Ar |Q|<Ar 

As in [31 [7] or [S], one can further deduce from p.3p the dissipation of a, b and c. In fact, let 
£1 > 0, £2 > be arbitrary constants to be chosen later. From (|3.3P f, and (H^Sja, it follows that 



'".9lP- 



d 
It 



V / V^a"c-A(9"{I-P}u)dx + A V ||V:,(9" 



lQ|<Af-l 



|Q|<Ar-i 



|a|<Ar-l ^ Vl"!^^ |Q|<Af-l 



(IX311 a and ([331)2 imply 
3 
dt 



3 



A E liv.a"fef<62 E llv.a"Kc]f+ E liv.a"i?|p 

Q|<Ar-l \|q|<7V-1 |a|<7V-2 



-I E ii5"{i-p}t 



E \\--'"d'^'-^ 



|Q|<Ar-l 
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It holds from (1331)2 and (l33|i that 



|Q|<Af-l"^*'^ \a\<N 



1q|<JV-1 |q|<JV-1 



^.11^ 



Define 



\a\<N—l 

3 

(3.6) + y y / (a,a"fe, + 9,o"6,;-^,5yV,-a"6)e,j(a"{i-p}u)da 



\a\<N-lij = l ' 

Ki V / V^d'^a-d'^bdx, 

|a|<N-l-^K' 



for some small constant ki > 0. Therefore, by taking ki > small enough and then letting ei > 
and £2 > be small enough, the dissipation of a, h and c can be obtained by the following inequality 

(3.7) l4"^i([/(i))+A Y. l|V.a"[a,6,c]f + ||af 

|a|<7V-l 

<^2 Y ||v.a"£;|p + ^( 5] ||a"{i-PMp+ Y W^-'^d^af 

|a|<iV-2 ^ \|a|<A' |a|<Ar-l 

where £2 > is still left to be chosen later on. 

The key part is to estimate the dissipation of E and B. We claim that 

(^•^) ^T Y. I d''E-d°'bdx + X Y ll^"^ll^ 

l<|a|<Ar-l "^^ 1<|q|<7V-1 

<£3 Y lia"V, Xi?||2 + ^ Y l|V.9"6f 

l<|a|<Ar-2 "^ |Q|<Ar-l 

+ C Y (l|V.9"[a,c]||2 + ||V,9"{I-P}"!P), 

l<|a|<Af-l 

with £3 > to be chosen, and 

(3.9) -4 Y f d^V^xB-d^Edx + X y Wd'^V^xBf 

dt ^-^ 03 ^-^ 

l<|a|<Af-2'"* 1<|q|<7V-2 

< Y w^x-^Ef+c Y ii^"^iP- 

l<|a|<Ar-2 i<\a\<N -2 

For this time, suppose that p.Sp and p.9p hold true. Define 

(3.10) f]^"'^(t/(i)) = - Y I d°'E-d°'bdx-K2 V [ d^^V^ x B ■ d^Edx 



l<|al<Ar-l ^ l<|Q|<Af-2' 
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for some constant K2 > 0. Then, by taking K2 > and further €3 > both small enough, it follows 
from KB and (Ejl) that 



(3.11) |4"'2([/(i)) + A Y. ii5"^ii' + A E ii^"^ii 

l<|a|<7V-l 2<|a|<Af-l 



<C||V,6f + C Y. (l|V.a"[a,6,c]||2 + ||V,9"{I-P} 

1<|q|<W-1 



u\\'] 



which is the desired dissipation estimate of the electromagnetic field E and B. Now, define 

(3.12) £:]^"(C/(i)) = E (liaNP + lia"[i?,i?]|p) + «4(4"'^(C/(t)) + «34"''(t^W)) 

|0!|<Af 

with constants K3 > and K4 > 0, where i5^"' {U{t)), £^' {U{t)) are defined in (|3.6[) and p.lOp . 
In the same way as before, by taking properly small constants K3 > 0, 62 > and K4 > in turn, 
p.gp follows from the linear combination of (|3.5p . (|3.7p and (|3.1ip . Moreover, it is easy to verify 
that £fj\U{t)) is the desired L^ energy functional satisfying PTT)) and 'D^}^{U{t)) is given by (fTS)) . 
Here, one has to check 

(3.13) E l|5"Sf<C E W^.y^Bf. 

2<|q|<W-1 l<|a|<Af-2 

In fact, by taking a with 2 < |q;| < A^ — 1 and using B = — A^r^V^; x Vj^ x i? due to V^; • i? = 0, 
it holds 

9"B = -a^A-^v^ X V,. X S = -a"^'''c),A-iV^ X V^ X B, 

for some 7i (1 < i < 3) with |7i| = 1. Since di/S.~^dj for any 1 < i,j < 3 is a bounded operator 
from LP to itself with 1 < p < 00, 

lia"B|i<c|i9"-T'V, xs|i. 

Hence, p.l3p follows from taking summation of the above inequality over 2 < |q;| < A^ — 1. 

Now, the rest is to prove ((3?8)) and ((3^ . Take a with 1 < |a| < iV - 1. By using ((3^ 9 to 
replace _E and then using p.4p i to replace dtE, one can compute 

(3.14) Wd'^Ef ^ [ d^E-d^Edx 

= [ d"'E-d"'[dtb + \/^{a + 2c) + \7^e{{I-P}u)]dx 
= ^ I d^E ■ d^bdx ~ I d'^dfE ■ d^bdx 



+ / d°'E-d°'[\7^{a + 2c) + V^e{{I-P}u)]dx 
-^ / d°'E-d°'bdx+ [ d^ib-V^x B)-d°'bdx 

dt Jr3 Jr3 

+ / d°'E-d°'[V^{a + 2c) + V^e{{I-P}u)]dx. 

JM3 



Then, ()3.8p follows from the above identity after taking summation over 1 < |a| < A^ — 1 and 
further applying the Cauchy-Schwarz inequality and integration by parts. In fact, it suffices to 
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consider the second term on the r.h.s. of p.l4|) . It can be estimated by 
Y^ f d°'{b~V^xB)-d"bdx= ^ ||5"6| 

\^.\^AT T"^ T ^ \ ^ \ ^ AT 1 



^ai,\\'^ 



l<|a|<Ar-l" 



l<|a|<JV-l 

Y^ I 5"V^ X B ■ d'^bdx + Y I d°'-"''\7^ x B ■ d°'+"''bdx 



l<|a|<iV-2 



\a\ = N~l 



c 



l<|al<Ar-2 "^ |a|<7V-l 

where as before < 63 < 1 is small to be chosen, and 7^ denotes a multi-index with |7,;| = 1 for 
some 1 < i < 3. To prove ([3^ . take a with 1 < |a| < TV - 2. By using ([5^ i to replace V^^ B 
and then using p.4p 9 to replace dtB, one has 

|ia"V:r X B||2 = ( d^V^ X B ■ d"V^ X Bdx = / a"V^ X B ■ d°'{dtE + b)dx 

Jr3 Jm.^ 



= — / a"V, X B ■ d°'Edx 
dt 



a"V, X dfB ■ d°'Edx - 



a"V, X B ■ d'^bdx 



a"V^ X B ■ d^Edx + / 9"V^ X V^ X £; • d^Edx + / 5"V^ x B • 9"6da;. 

Then, p.9p follows by applying integration by part to the right-hand second term of the above 
identity, using the Cauchy-Schwarz inequality and then taking summation over 1 < |a| < N — 2. 
The proof of Theorem 1 1.1 1 is complete. D 

3.2. L^ time-frequency functional and its optimal dissipation rate. In this subsection we 
shall prove Theorem 11.21 as well as Corollarv ll.il and Corollarv ll.2l For that, wc need to consider 



the solution U = [u,E,B] to the Cauchy problem p.ip - ()3.2p in the Fourier space R|. By taking 



the Fourier transform in x from (|3.ip 1 , (|3.3p and p.4p , one has 



(3.15) 



(3.16) 



and 



(3.17) 



dtu + i^-ku~ ^Mi/2 .E^Lu + h, 

dto, + ik ■ b = 0, 

dtb + ik{a + 2c) + ik ■ e({I - P}u) - £; = 0, 

dtc +-ik-b+ -ik ■ A({I - P}u) = 0, 



dtQij ({I - P}u) + ikibj + ikjbi ~ -S^jik ■ b 
^ a* A, ({I - P}u) + ihc = A.,ie + h), 
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^s,jik-A{{i-p}u) = e,jie + h), 



dtE - ik X B ^ -6, 

dtB + ikx E = 0, 

ik ■ E = a, ik ■ B = 0, 



where i is given by 

i = -ik ■ ^{I - P}m -I- Lu. 
These equations above are ones to be used through this subsection. 
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Proof of Theorem II. 2t It is similar to the proof of Theorem 11.11 The difference is that all 
calculations are made in the Fourier space. Thus, some details in the following proof will be 
omitted for simplicity. First of all, as in [3], on one hand, from (|3.15p and p.l7p . one has 

(3.18) \dt{\\u\\l. + \[E, B]\^) + Ai|^i/2|i _ p}^||2^ < C||^-i/2g|ii., 
and on the other hand, from p.l6p . the following three estimates hold true: 

(3.19) dt'Kiikc \ A({I - P}u) + X\k\'^\c\'^ 

<e,\\k- hr + -(1 + IfcHlllI - V}u\\l. + C\\u-^'^g\\l., 



(3.20) dt^^{ihhj+ikjk - i;Sijik ■ b \ e.y({I- P}m) + X\k\'^\b\'^ 

< e,\kna, c]|2 + e,-^\E\^ + -(1 + \k\')\\{I - P}u\\l. + C\W-'/'g\\l., 

i + \K\^ £2 5 5 

and 

(3.21) dt^niika \ b) + A(l + \k\'^)\a\'^ < \k ■ b\^ + C\k\^\c\^ + C\k\'^\\{I - P}u\\l2 



^V 



where constants < ei, £2 < 1 are to be chosen. Define 
(3.22) £''-'\lJ{t)) = ^±-^{(^kc \ A({I - P}^) 



y {ikibj + ikjbi — —Sijik ■ b \ O^ ({I — P}m) + Ki{ika \ b)} 



3 

for a constant ki > 0. One can take ki and then ei both small enough such that the sum of p.l9p . 
(P:^ and Ki X ((X^ gives 

(3.23) d,£''--\U{t)) + A J^|[a,S,c]|2 + \af 

For estimates on the dissipation of E and B, it is straightforward to deduce from p. 15^ 9 and 
p.l6p the following two identities 

(3.24) - dt{k X E\kxb) + \kx E\'^ = \kxb\^ -{ikxkx B\kxb) 

+ {kxE\ikx{k- e({I - P}u))), 
and 

(3.25) - dtiik X B \ E) + \k x B\^ = |fc x E\^ + {ik x B \ b). 

By applying the Cauchy-Schwarz to p.25p and then multiplying itby|fcp/(l + |fcp)^, one has 

-^'T^TW¥^ X B I i?) + A^^-^^^ < (i + |fc|2)3 + ^(l + |fcp)3' 
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which imphes 

Similarly, after dividing p.24p by (1 + |fcp)^ and then using Cauchy-Schwarz, 
(3 27) dMk xE\kxb) ^ ^^ \k X E\^ 



(l + |fc|2)2 (l + |fc|2)2 

|A;x6|2 |fcxfcx^|-|fcx6| |fc|^|9({I ~ P}m)P 

- (l + |fcP)2 + (l + |fc|2)2 + (l+|fc|2)2 

^ |fc|2|fcx^|2 C |fc|2 

-'' (i + |fc|2)3 +^TT^I^I +C11II-PHL., 



where we used the inequality 



|fcx A:xB| • |fcx 6| \k x k x B\^ C \k x b\^ 

< £3- ' 



(l + |fc|2)2 - ^ (1 + |A:|2)3 e3l + |fc|2 

for an arbitrary constant < 63 < 1. Then, in terms of p.26p and p.27p . let us define 
.o^n^ ^lin,2rm.^^ ^{kxE\kxb) \k\^D\{ik xB\E) 

(3-28) £ {U{t))- (l+|fc|2)2 -«2 (i+|fc|2)3 ' 

where K2 > is chosen small enough such that 



|fcxig|2 |fc|4|i?|2 ^ |fc|2|6|2 

'(l + |fc|2)2^ (1 + |A:|2)3- i+|fc|2'-!ii- -;-iiL? 



(3.29) a,£'-^2(^(^)) ^ A '",,,' , + A-P-^:^ < C^^^^ + C!1{I - P}h||2, 



Here, we used \k x B\ = \k\ • |fc x i?| due to fc • i? = 0. 

Now, in terms of ((XT^ . ([5:^ with < 63 < 1 and (jX^ . we define 

(3.30) £""(;7(t)) = llTilli. + |[^, Bf + K4 (£:""-^(?7(i)) + K,8'"'^\U{t))) , 

where f""'! ([>(<)), f""'2((7(^)) ^re denoted by ([X^ and ([X^ . and K3 > 0, 62 > and K4 > 
are chosen in turn small enough such that (jl.lOp and (|1.12p hold true and I?''"(f7(t)) is given by 
(jl.lip . Here, notice that we used 

|A:|2|^|2 = \k ■ E\^ + \kx E\^ =^ |a|2 + \k x E\^. 

The proof of Theorem 11.21 is complete. D 

Proof of Corollary 11.11 and Corollary II. 2t First, (J1.13p in Corollary 11.11 immediately results 
from (|1.12p by noticing 

ViU{t,k))>X^^^^SiU{t,k)) 
due to the definitions ([TTU)) . ([rTTj) oi£{U{t, fc)), V{U{t, k)). For Corollary [TT^ it suffices to define 

(3.31) £'^\U{t,k))^ Y. n^it,k))= ^ £{{tkrU{t,k)). 

\a\—m |a|— m 

Since the system p.ip satisfied by f/ = [u, E, B] is linear, it is easy to see from Theorem 11.21 that 
£ll^{U{t,k)) satisfies ([TTi)) and p:TB|) with V^^'^{U{t,k)) given by dUg). This hence completes 
the proof of Corollary 11.21 D 
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3.3. Time-decay estimates. In this subsection we turn to the proof of Theorem 11.31 Theorem 
II. 31 actually follows from (|1.13p in Corollarv ll.il However, we would rather provide a similar much 
more general result. Recall the definition of A(t) as in p.l8|) . In what follows, for the Fourier 
transform U{t, k) oi U = [u, E, B], we set 

(3.32) \U{t,k)\^\\u\\L2 + \[E,B]\ 

for simplicity. 

Lemma 3.1. Assume that for any initial data Uo, the linear homogeneous solution U{t) ~ A{t)Uo 
obeys the pointwise estimate 

(3.33) |C>(i,fc)| <Ce-'^('=)*|f7o(fc)| 

for all t > 0, k ^ M.'^ , where (f){k) is a strictly positive, continuous and real-valued function over 
k G M.^ and satisfies 

rO(l)|fc|'^+ as|fc|^0, 

(3.34) </>(fc) ^ <^ 

iO(l)|fc|-'"- as|fc|^oo, 

for two constants cr_ > cr+ > 0. Let m > be an integer, 1 < p,r < 2 < q < oo and cr > 0. Then, 
U{t) = A(t)[/o obeys the time-decay estimate 



3_/ 1_1 



\ n Q I 



v--^ii,-^™+['^+3(i 



(3.35) ||V;"[/(i)||2, <C(l+i) "+^- -' '^\\V^\\z,+C{\+t) -IJV, 
for any t>0, where [■]-^- is defined in (jl.22p . 

Proof. Take a constant R > 0. From the assumptions on (j){k), it is easy to see 

A|fc|'^+ if |fc| < R, 

X\k\-''- if |fc| > R. 
Take 2 < q < oo and an integer ?7i > 0. From Hausdorff- Young inequality. 



Ui 



0\\Zr 



0(fc)> 



(3.36) ||V™C/(f)b, <C 



z„, 



<C 



\k\me-m'+tfj^ 



c 



m „ — A I A: I 



\kre 



Z,,(\k\<R) 

where ^ + ^ = 1 . For /i , by the definition of the norm 1 1 • \\z ,, 



'Uo 



z,'m>R) 



h + h, 



h^C 



|^|mg-A|fcr+t^^ 



Ll(Lo'm<R)) 



c 



|fcre-^l'=l'^*[^o,So] 



Li'(\k\<R) 



Here, note that since 1 < 9' < 2, from the Minkowski inequality. 



IfcPe- 



-Aifcr+t- 



MO 



< 



LliLo' {\k\<R)) 



\kre 



-Ai/tr+t - 



Wo 



Lt' (\k\<R-Ll) 



va-\\k\''+t\ 



Hence, we arrive at 



h<C 



\kre 



-Alfer+ir 



uo||l?,-E'o,So 



\kre 



Li' (\k\<R) 



uoWlI 



Li' {\k\<R) 



Take 1 < p < 2. Further using the Holder inequality for 4- — ^ , ? + -h with p' given by i + -^ = 1, 



h<C 



\kre 



-Aifcr+t 



V 1 

Lp'-i' (|fc|<-R) 



\ua\\L2,Eo,Bo] 



Lp'{\k\<R) 
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Here, the right-hand first term can be estimated in a standard way |12j as 

>l7icri 3 /1 1 \ m 



Lp'-i' (|fe|<-R) 
1 

by using change of variable kt''+ -^ k, and the right-hand second term is estimated by Minkowski 
and Hausdorff- Young inequahties as 



|wo||l;,£^o,So] 



— !l"0|lLP'(|fe|<H;L2) + 

< \\^o\\Ll{Lp'{\k\<R)) + 



Eo,Bf) 



Lp'{\k\<R) 



Lp'{\k\<R) 

< C(||lio||L^(LS) + \\[Eo,Bo]\\lp) = C\\Uo\\z,, 
where the Minkowski inequality was validly used due to p' > 2. Therefore, for /i, one has 

h<c{i + ty^^-'--'>'^\\Uo\\z,. 

To estimate I2, take a constant cr > so that 



h = C 



|fc|"e 



-Xlkl^'-t 



Un 



< C sup Ifcr'^e-^''^'"""* Ikr+'^Uo 

Z,'i\k\>R) \k\>R 



Z„'i\k\>R) 



Here, the right-hand first term decays in time as 



sup |fc|-''e-^l'^'l " * <C{l + ty 

\k\>R^ 



We estimate the right-hand second term as follows. Take 1 < r < 2 with - + ^ = 1 and take a 
constant e > small enough. Then, similarly as before, from Minkowski and Holder inequalities 
for J7 = '' r? + ^, one has 

q' r q r ' 



< 



< 



Z,,{\k\>R) 
|^|-3(1+.)^ 



\kr+''[\\uo\\LlEo,Bo] 



Li'{\k\>R.) 



\k\' 



i+o-+3(l+e)^-74- 



\i'-o\\L^,EQ,Bo] 



L-'{\k\>R.) 



<a 



\k\ 



'n-l-k+3(i-i)]+r 



luollLhEo^Bo] 



L'-'(|fc|>_R) 



Here, by Minkowski inequality due to q' > 2 once again and further by Hausdorff- Young inequality, 

,m+[£T+3(i-i)]+r 



< 



\k\ 



™+k+3(i-i)]+ 



Uq 



LjiL'-' {\k\>R.)) 



\k[ 



[En, Bo] 



L^'{\k\>R} 



<ciivr'^"^^^-'^'-c/oib. 



Then, it follows that 



\kr+^uo 



2„'(|fc|>-R) 



<ciivr"^+'^^~'^'^c/oiu.. 



Thus, I2 is estimated by 



7.<c(i+t)--iivr'^+^'-^''-f/oiiz.. 



Now, p.35p follows by plugging the estimates of /i and I2 into p.36p . This completes the proof 
of Lemma 13.11 D 
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Proof of Theorem [mi To prove pr20| . by letting /i = and using Corollary [LTI 



A|fc|'» 



£{U{t,k))<e (^+i''i-'y' £{Uo{k)) 

for any t > and fc S M'^, where we have set U = U^ for simplicity. Due to (|1.10p and p.32p . 
£{U{t,k)) - |C/(t,fc)P holds so that 

\U{t,k)\ < Ce"^aWF*|f7o(fc)| 

for any t > and fc G M'^ . This shows that corresponding to (|3.33p and (|3.34p of Lemma 13. 1[ one 
has the special situation 

'^^') = 2(l + |fcP)3 

with (7+ = 4, cr_ = 2. Thus, one can apply Lemma [3. II to obtain p.20p from p.35p . 

To prove (|1.2ip . we let Uo ~ 0, and also set U = U^' for simplicity. Note that (|1.13p implies 



\U{t,k)\^<C f e ^^oTTW^* "\\iy-^^^h{s)\\l2ds 
Jo '' 

for any i > and /c G M'^. One can again apply Lemma [3. II with g = r = 2 so that p.2ip follows. 
The proof of Theorem 11.31 is complete. D 

4. Nonlinear system 

In this section we are concerned with the global existence of solutions to the Cauchy problem 
p.3p - p.4p of the reformulated nonlinear Vlasov-Maxwell-Boltzmann system. We first devote our- 
selves to the proof of some uniform-in-time a priori estimates on the solution. In what follows, 
U = [u, E, B] is supposed to be smooth in all arguments and satisfy the system (|1.3I) over < t <T 
for some < T < oo. 

Lemma 4.1. Under the assumption that supg^j^-^ || (a + 2c)(i)||Loo is small enough, there is 
£N{U{t)) satisfying (|1.23|) such that 

(4.1) ^£N{U{t)) + \VN{U{t)) < C{£N{U{t)f'^ + £N{U{t)))VN{U{t)) 
at 

for any < t < T, where I?Ar([/(t)) is defined in (|1.24p . 
Proof. First, the zero-order energy estimate implies 

(4.2) ~i\Wr + II [E, B]r- j \b\\a + 2c)dx + \\\v^'^{l V}uf 

< C{£N{U{t)f' + £N{U{t)))VN{U{t)). 

Here and hereafter, when £N(U(t)) occurs in the right-hand terms of inequalities, it means an 
equivalent energy functional satisfying (|1.23p and its explicit representation will be determined 
later on. In fact, from the system (|1.3p . 

(4.3) ~{\\uf + \\[E,B]f) + X\\iy^^^{I-P}u\\^< j j uT{u,u)dxd£, + ]- jh-Eu^dxdS,. 

Here, for the first term on the r.h.s. of (|4.3p . it is a standard fact as in [9] or [10] that it is bounded 
by C£NiU{t)y^'^VN(U(t)). Since all estimates on terms involving the nonlinear term T{u,u) in 
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the following can be handled in the similar way, we shall omit the details of their proof for brevity. 
The right-hand second term of (|4.3p can be estimated as in [71 Lemma 4.4]: 

- ff ^- Eu^dxd£. ^- jh- E\Vu\^dxd^ ^ [[ ^' -^P"{I - P}udxd^ 

+ \jj i-E\{l-V}u\Hxdi. 
Here, it is easy to see 

^ • EI>u{l - I>}udxd£, + i jh- E\{1 - 'P}u\^dxdS, 

<C\\E\\H^{\\V,[a,b,c]f + \\v^'^{l-V}uf)<C£N{U{t)Y'^VN{U{t)). 

And, as in [7J Lemma 4.4], one can compute 

\ II ^- E\Pu\^dxdS, = I E-b{a + 2c)dx, 

where by replacing E by equation (J2.9p 9. using integration by part in t and then replacing dt{a + 2c) 
by equations p. 9^ 1 and ()2.9p -;. gives 

E ■b{a + 2c)dx = 7:-r I l^l^(" + 2c)da; 

+ / l^n^V, ■ & + ^V, • A({I - V}u) - ^E ■ b]dx 
+ / [^xia + 2c) + V:r6({I - P}u)] • b{a + 2c)dx 
Ea ■ b{a + 2c)dx - b x B ■ b{a + 2c)dx. 



Note b X B ■ b = 0. Hence, it follows that 

^JJ^.E\Pu\'dxd^<~J\b\\a + 2c)dx 

+ C\\[a,b,c,E]\\Hi{\\^Aa,b,c]\\^ + \\WAl-P}uf) + C\\E\\\\WM\\\^x[a,b,c]f 

-\jtj '^''^" + ^""^'^'^ + C{£N{U{t)Y'^ + £N{U{t)))VN{U{t)). 

Collecting the above estimates and putting them into ()4.3p proves (|4.2p . 

Next, for the estimates on all derivatives including the pure spatial derivatives and space- velocity 
mixed derivatives, one has 



l<|Q|<Af ^<\a\<N 



U\\^ 



< C£N{U{t)f'^VN{U{t)), 
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and 



(4-5) ^1E^^ E II5|{i-phP + a ^ ||.i/2a^ni_p},||2 

k=l \l3\=k |^|>1 

|a| + |0|<Af |a| + l/3|<W 

Q|<Af-l \a\<N 



where Cfc (1 < fc < A'') are strictly positive constants. Since the nonlinear term g takes the form as 
in (|2.8p . the proof of (|4.4p and (|4.5p is almost the same as in [^ for the case of the Vlasov-Poisson- 
Boltzmann system and thus details are omitted for brevity. 

Finally, the key step is to obtain the macroscopic dissipation of a, b, c and E, B for the nonlinear 
system (|1.3p . This is similar to the proof of Theorem 11.11 in Subsection 13.11 for the linearized 
system. Here, the additional efforts should be made to take care of all quadratically nonlinear 
terms in g defined by (|2.8p . But, these nonlinear estimates once again are almost the same as in 
[7] for the case of the Vlasov-Poisson-Boltzmann system so we omit details for brevity. Thus, we 
have the following estimates. Recall the definitions p.6p and p.lO[) of two interactive functional 
£^' {U{t)) and £™' {U{t)), where ki > 0, K2 > in p.6|) and ()3.10p are sufficiently small. It 
turns out that 



(4.6) j^[£p\U{t)) + ^s£f-\U{t)))+xY. \\d"af + X ^ \\d"[b,c]f 

\a\<N l<|a|<Ar 

+ A Y. Il^^-Bf + A Y. ll^"^f 

l<\a\<N~l 2<|a|<W-l 

<C Y W{I-P}uf + C£N{U{t))VN{U{t)), 

\a\<N 



where K3 > is small enough. This is the desired estimate on the macroscopic dissipation. 
Now, we are in a position to prove (|4.1I) . Let K4 > be taken as in (|3.12p . Define 



N 

(4.7) SNiUit)) = f]^"(C/(i)) - / \b\\a + 2c)dx + ^5 E ^^ E H^^i^ " P^^U'' 

■^ k=l \f3\=k 

\a\ + \l3\<N 
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where K5 > is a constant to be chosen. Here, notice that by recalhng the definition (|3.12|) of 
^N^iUit)), £N{U{t)) can be rewritten as 

(4.8) £N{U{t)) 

= ^ {\\d"uf + W[E,B]f)- f\b\\a + 2c)dx 

\a\<N ■' 

+K4< y / V^d°'c-K{d°'{l-V}u)dx 

+ y. f (^^^"^3 + ^j^^b, - h,jV, ■ d"b)e^jid"{i - p}u)dx 

\a\<N~l 

+Ki V / V^d^a-d^bdx} 

J - y f d°'E ■ d°'bdx -K2 Yl I ^"' 

AT 

+K,Y,cu Y. ii5|{i-PMp. 

fc=l \l3\=k 

\a\ + \l3\<N 

Due to smallness of supo<f<T II (a + 2c)(i)||L= by the assumption, from ()4.7p and ()1.7p . it is easy 
to see 

N 

EN{U{t))^El^{U{t)) + YCk J2 l|5^{I-P}«f, 
which further imphes 



+K4K3 { y / d^E ■ d^bdx -K2 y / 5"V:, X B ■ d^Edx 



k=l \p\=k 

|a| + l;3|<Af 



AT 

2 



U\ 



8M{u{t)) ^ i|u(i)iii.(^«) + ii[ii;(t),i3(i)]|l?,« +E^^ E ll^^"{i-P} 

fe=l \0\=k 

\<x\ + \l3\<N 

^||u(i)||^„ +\\[E[t),B{t)]\\l.. 

Thus, £]\[{U{t)) satisfies (|1.23p for any K5 > 0. Moreover, by taking K5 > small enough, the 
summation of (g^l), (|0)) . KziX (|TB)) and then K^x ff^ leads to (|ITT|) with 2?Ar(t/(i)) defined in 
(jl.24p . This completes the proof of Lemma HTTl D 

Proof of Theorem II. 4t Let us consider the uniform-in-time a priori estimates of solutions under 
the smallness assumption that 



sup 

0<t<T 



[\\umi^^^ + \\[Eit),Bimi.)<s 



for a sufficiently small constant S > 0. This smallness assumption also implies that supo<t<'r II (a + 
2c)(i)||L=^ is small enough since A^ > 4. Thus, it follows from Lemma ITT] that there are £N{U{t)), 
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VN{U{t)) defined in (gl]) and (fLM]) such that (fr23| holds true for EN{U{t)) and 

^EMiUit)) + AI?^([/(i)) < C{5''^ + <5)I?^(C/(i)) 
at 

for any < t < T, that is 

£N{U{i)) + X f VN{U{.s))ds < £n{Uo) 
Jo 

for any < t < T, since (5 > is small enough. Now, the rest proof follows from the standard 

process by combining the above uniform-in-time a priori estimates with the local existence as well 

as the continuity argument as in [TT] or [7] under the assumption that En{Uq) is sufficiently small, 

and details are omitted for simplicity. The proof of Theorem ll.4l is complete. D 

We conclude this paper with a discussion about the large-time behavior of solutions to the non- 
linear system. Although Theorem 11.41 shows the global existence of close-to-equilibrium solutions 
to the Cauchy problem (|1.3p - (|1.4p of the nonlinear Vlasov-Maxwell-Boltzmann system, the decay 
rate of the obtained solution remains open. This issue has been studied in [5] for the case of 
two-species. However, the approach of [S] by applying the linear decay property together with the 
Duhamel's principle to the nonlinear system can not be applied to the case of one-species here. 
Let us explain a little the key difficulty in a formal way. In fact, the linear system in one-species 
decays as (1 -I- t)~^^^ which is slower than (1 -I- t)~^/^ in two-species as pointed out in Table 1. 
Thus, in one-species case, the quadratic nonlinear source decays as at most (1 -I- i)~^/^, and if the 
Duhamel's principle was used, the time-integral term generated from the nonhomogeneous source 
decays as 

rt 

(l + i-s)"S(i + s)-Sds < C{l+t)-i. 

10 

Thus, the bootstrap argument breaks down and one can not expect the solution to decay as 
(1 -|- t)""^/^ in the nonlinear case. Therefore, the study of the large-time behavior for the one- 
species Vlasov-Maxwell-Boltzmann system becomes much more difficult than for the two-species 
case as in [5]. 
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